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ABSTRACT: 

This  report  summarizes  the  results  corresponding  to  a  series  of  dynamic  stress  analyses  performed  on  an 
overflow  and  a  nonoverflow'  monolith  of  Folsom  Dam,  California.  The  monoliths  selected  for  analysis 
were  monoliths  14  and  21,  and  Aey  were  chosen  as  representative  of  critical  overflow  and  nonoverflow 
configurations.  The  analyses  were  conducted  using  the  computer  program  EAGD-84,  developed  by  Pro¬ 
fessor  A.  K.  Chopra  and  co-workere  at  the  University  of  California  at  Berkeley.  The  analyses  accounted 
for  dam-foundation  interaction  and  water  compressibility  effects,  and  considered  two  lo^ng  scenarios: 
(a)  horizontel  component  of  the  ground  motion  and  (b)  horizontal  and  vertical  components  of  the  ground 
motion.  The  results  are  reported  in  terms  of  the  maximum  expected  stresses  under  extreme  loading  condi¬ 
tions  (maximum  credible  earthquake,  MCE),  and  they  are  compared  wdth  those  results  obtained  from 
analyses  performed  by  engineers  at  the  Sacramento  District  using  a  response-spectrum-based  procedure 
also  developed  bj^  Professor  Chopra.  The  influence  of  the  vertical  component  of  the  ground  motion  on  the 
magnitude  and  distribution  of  the  d^mamic  stresses  is  investigated.  The  levels  of  tensile  demands  pre¬ 
dicted  by  this  study  indicate  that  localized  damage  in  the  form  of  tensile  cracking  is  likely  to  take  place  in 
limited  areas  of  monoliths  14  and  21  when  subjected  to  extreme  seismic  loading,  but  this  is  still  accept¬ 
able  seismic  performance  under  MCE  conditions. 
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Conversion  Factors, 
Non-SI  to  Si  Units  of 
Measurement 


Non-SI  units  of  measurements  used  in  this  report  can  be  converted  to  SI  imit 
as  follows: 


Multiply 

By 

To  Obtain 

cubic  feet 

0.028317 

cubic  meters 

Fahrenheit  degrees 

(F-32)/1.8 

Celsius  degrees 

feet 

0.304800 

meters 

gallons 

0.00378 

cubic  meters 

grams 

0.001 

kilograms 

inches 

25.4 

millimeters 

miles 

1.609 

kilometers 

ounces 

0.00002957 

cubic  meters 

pounds  (force)  per  square  inch 

0.006894757 

megapascals 

pounds  (mass)  per  cubic  foot 

16.01846 

kilograms  per  cubic  meter 

square  inches 

6.4516 

square  centimeters 

Preface 


This  report  describes  a  research  study  developed  by  the  U.S.  Army  Engineer 
Research  and  Development  Center  (ERDC)  consisting  of  dynamic  stress  analyses 
performed  on  monoliths  14  and  21  of  Folsom  Dam,  California.  These  monoliths 
were  chosen  as  representative  of  critical  overflow  and  nonoverflow  configura¬ 
tions.  The  analyses  were  conducted  using  the  computer  program  EAGD-84, 
developed  by  Professor  A.  K.  Chopra  and  co-workers  at  the  University  of  Cali¬ 
fornia  at  Berkeley.  The  results  from  these  analyses  complement  the  information 
provided  by  other  previous  and  current  technical  studies,  and  it  is  expected  that 
they  will  contribute  toward  a  more  accurate  assessment  of  the  seismic  perfor¬ 
mance  of  this  critical  structure.  The  study  was  part  of  the  Civil  Works  Direct- 
Allotted  354  Infrastmcture  Technologies  R&D  Program,  Work  Unit  33013, 
“Rehabilitation  of  Concrete  Hydraulic  Structures  Under  Extreme  Loads,”  which 
is  sponsored  by  Headquarters,  U.S.  Army  Corps  of  Engineers. 

This  publication  was  prepared  by  personnel  from  the  ERDC,  Geotechnical 
and  Structures  Laboratory  (GSL),  Vicksburg,  MS.  The  research  described  herein 
was  conducted  by  Dr.  Enrique  E.  Matheu  and  Mr.  Raju  V.  Kala,  Geotechnical 
and  Earthquake  Engineering  Branch  (GEEB),  and  Dr.  Robert  L.  Hall,  Chief, 
Geosciences  and  Structures  Division  (GSD),  GSL.  Dr.  Matheu  prepared  this 
publieation  under  the  general  supervision  of  Dr.  David  W.  Pittman,  Acting 
Director,  GSL;  Dr.  Robert  L.  Hall,  Chief,  Geosciences  and  Structures  Division, 
GSL;  and  Dr.  Joseph  P.  Koester,  Chief,  GEEB,  GSL. 

At  the  time  of  publication  of  this  report,  COL  James  R.  Rowan,  EN,  was 
Commander  and  Executive  Director  of  ERDC,  and  Dr.  James  R.  Houston  was 
Director. 
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Introduction 


Objective 

The  objective  of  the  technical  work  documented  in  this  report  is  to  evaluate 
the  dynamic  response  of  Folsom  Dam  in  the  event  of  the  Maximum  Credible 
Earthquake  (MCE).  This  technical  study  included  tw'o-dimensional  dynamic 
stress  analyses  of  representative  sections  of  the  concrete  dam.  The  U.S.  Army 
Engineer  District,  Sacramento  (SPK),  commissioned  this  study  to  determine  the 
expected  stmctural  performance  of  the  main  concrete  dam  under  extreme  loading 
conditions.  An  important  objective  is  to  compare  the  results  obtained  in  this 
study  by  conducting  time-history  finite-element  analyses  of  two-dimensional 
sections  with  the  corresponding  results  calculated  by  the  Sacramento  District 
using  response-spectrum-based  analyses  of  three-dimensional  models.  The 
results  and  conclusions  presented  in  this  report  are  intended  to  complement  the 
information  provided  by  other  previous  and  current  technical  studies,  and  it  is 
expected  that  they  will  contribute  towards  a  more  accurate  assessment  of  the 
seismic  performance  of  this  critical  structure. 


Scope 

This  report  presents  the  results  of  a  series  of  finite -element-based  dynamic 
stress  analyses  of  selected  two-dimensional  sections  of  Folsom  Dam.  The 
selected  sections  corresponded  to  monoliths  14  and  21  of  the  main  concrete  dam. 
Sacramento  District’s  technical  persoimel  selected  these  two  monoliths  (14  and 
21)  as  representative  of  the  most  relevant  dynamic  response  characteristics  of 
critical  overflow  and  nonoverflow  sections  of  the  concrete  dam,  respectively.  The 
analyses  are  performed  using  well-recognized  linear  time -history  response 
procedures  and  the  results  include  the  maximum  expected  stresses  under  extreme 
loading  conditions  (MCE),  accounting  for  dam-foundation  interaction  and  water 
compressibility  effects.  The  results  from  this  numerical  study  complement 
previous  analyses  performed  by  the  Sacramento  District  using  the  simplified 
analysis  procedure  developed  by  Professor  Anil  K.  Chopra  (Chopra  1978,  Fenves 
and  Chopra  1986,  Chopra  and  Hanchen  1989). 
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Background 

Project  characteristics 

Folsom  Dam  is  located  on  the  American  River,  about  20  miles  northeast  of 
the  city  of  Sacramento,  CA.  The  reservoir  is  used  not  only  to  provide  flood  con¬ 
trol,  but  it  also  serves  for  iirigation  and  power  generation  purposes.  The 
maximum  height  of  the  gravity  dam  section  is  340  ft'  with  a  crest  length  of  about 
1,400  ft.  It  consists  of  28  monoliths,  50  ft  wide  each,  and  its  construction  was 
completed  in  1956.  The  monoliths  were  constructed  in  5-ft  lifts  and  are  founded 
in  hard  granodiorite  rock.  A  2-  to  10-ft-thick  shell  of  concrete  with  relatively 
high  cement  content  (rich  concrete)  was  placed  along  the  upstream  and 
downstream  faces  of  the  monoliths  from  the  base  to  the  crest.  A  concrete  mix 
w'ith  lower  cement  content  (lean  concrete)  was  placed  throughout  the  rest  of  the 
dam  section  (Hall  et  al.  1989). 

Ihe  motivation  of  the  Folsom  Dam  Outlet  Works  Modifications  Project  is  to 
reduce  the  risk  of  flooding  to  the  city'  of  Sacramento  posed  by  the  American 
River.  Specifically,  the  objective  of  this  project  is  to  increase  the  current  outflow' 
capability  from  26,000  to  1 15,000  ftVsec  at  the  reservoir  w'ater-surfece  elevation 
corresponding  to  the  spillw'ay  crest.  This  increased  release  capacity  w'ill  lead  to  a 
more  efficient  use  of  the  flood  reservation  space,  and  w'ill  increase  the  overall 
level  of  downstream  flood  protection  by  reducing  the  probability  of  flooding  in 
Sacramento  in  any  one  year  from  one  chance  in  100  to  one  chance  in  140.  The 
project  basically  consists  of  enlarging  the  four  existing  upper  tier  river  outlets 
from  5  ft  w'ide  by  9  ft  high  to  9.33  ft  wide  by  14  ft  high,  constructing  two  new' 
upper  tier  river  outlets  of  the  same  size,  and  enlarging  the  four  existing  lower  tier 
river  outlets  from  5  ft  wide  by  9  ft  high  to  9.33  ft  wide  by  12  ft  high  (W ong  et  al. 
2002). 


Previous  analyses 

The  Corps  conducted  a  series  of  comprehensive  seismic  evaluation  studies  of 
the  dam  in  the  late  1980s.  These  studies  included  stress  and  global  stability 
analyses  of  critical  sections  of  the  main  concrete  dam  (Hall  et  al.  1989).  The 
MCE  w'as  defined  as  an  event  of  magnitude  6.5  at  a  source -to-site  distance  of 
about  15  km,  on  the  eastern  branch  of  the  Bear  Mountains  fault  zone.  The 
corresponding  peak  ground  acceleration  (PGA)  for  the  horizontal  direction  w'as 
defined  as  0.35  g,  and  two  different  horizontal  ground  acceleration  time  histories 
w'ere  selected  for  the  dynamic  analyses.  The  corresponding  vertical  ground 
acceleration  time  histories  were  generated  by  increasing  the  frequency  content  by 
a  factor  of  1 .5  and  reducing  the  amplitudes  by  a  factor  of  0.60. 

The  dj'namic  stress  analyses  were  performed  using  the  computer  program 
EAGD-84.  Preliminaty  analyses  of  several  cross  sections  revealed  that  the  tallest 
nonoverflow'  monolith  represented  the  critical  section.  The  results  from 


'  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI  units  is  presented  on 
page  viii. 
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additional  analyses,  considering  different  values  of  foundation  modulus  (5.8,  7.9, 
and  1 1.0- 10^  psi)  and  wave  reflection  coefficient  (0.75,  0.79,  and  0.82)  as  well  as 
all  possible  sign  combinations  for  the  ground  motion  pairs,  indicated  that  the 
stresses  due  to  the  MCE  were  within  a  range  deemed  acceptable.  The  most 
demanding  condition,  characterized  by  the  largest  values  of  the  foundation 
modulus  and  the  wave  reflection  coefficient,  resulted  in  a  maximum  principal 
stress  that  reached  a  peak  value  of  871  psi.  This  peak  value  occurred  on  the 
downstream  face,  near  the  lower  end  of  the  circular  transition. 

This  study  also  included  a  global  stability  analysis  of  the  tallest  nonoverflow 
monolith.  This  analysis  was  performed  using  a  seismic  coefficient  of  0.15  to 
pseudo-statically  represent  the  lateral  inertial  forces  acting  on  the  monolith.  The 
shear  strength  values  used  in  the  analysis  were  based  on  typical  values  and 
included  a  significant  contribution  from  cohesion  at  the  dam-foundation 
interface.  The  computed  factors  of  safety  against  sliding  were  considered 
acceptable,  but  some  of  the  underlying  assumptions  would  not  be  presently 
considered  valid  based  on  current  guidelines  for  analysis. 

As  part  of  the  Folsom  Dam  Outlet  Works  Modification  Project,  the  seismic 
performance  of  Folsom  Dam  was  evaluated  by  global  stability  and  stress  analyses 
of  selected  monoliths  performed  by  the  Sacramento  District.  The  earthquake 
excitation  was  defined  according  to  probabilistic  and  deterministic  seismic 
hazard  analyses  based  on  the  most  current  data  available  (URS  Corporation 
2001;  URS  Corporation  2002).  Stress  analyses  of  monoliths  13,  14,  and  21 
consisted  of  the  application  of  Chopra’s  simplified  response  spectmm  procedure. 
This  approach  was  used  to  determine  two  sets  of  earthquake  lateral  forces 
associated  with  the  imposed  upstream-downstream  base  motion  (the  method  was 
later  expanded  to  consistently  include  the  effects  of  vertical  ground  motion).  The 
first  set  of  forces  corresponds  to  the  lateral  inertia  forces  associated  with  the 
fundamental  vibration  mode  of  the  two-dimensional  model  of  the  monolith, 
considering  dam-water  interaction,  water  compressibility,  reservoir  bottom 
absorption,  and  dam-foundation  rock  interaction  effects.  The  second  set  of  forces 
corresponds  to  the  lateral  earthquake  forces  associated  with  higher  vibration 
modes  of  the  system.  An  underlying  assumption  in  the  simplified  procedure  used 
to  determine  the  earthquake  forces  is  that  the  monolith  is  supported  on  a 
horizontal  surface  of  underlying  flexible  foimdation  rock.  This  is  certainly 
appropriate  for  monoliths  14  and  21,  which  exhibit  a  dam/foundation  contact 
with  approximately  uniform  elevations.  However,  in  the  case  of  monolith  13,  the 
dam/foundation  contact  varies  by  as  much  as  35  ft  in  elevation,  and  it  was 
necessary  to  assume  a  horizontal  base  based  on  the  average  monolith  height. 

The  earthquake-induced  lateral  forces  were  applied  to  three-dimensional 
(3D)  finite  element  models  of  the  spillway  monoliths,  developed  using  the 
computer  program  GT  STRUDL.  Taking  advantage  of  sjonmetry  conditions, 
detailed  3D  models  were  developed  for  only  half  (25  ft)  of  the  corresponding 
width,  and  they  included  the  corresponding  conduits  and  chamber  openings.  In 
the  analysis  of  the  spillway  monoliths,  two  different  models  were  generated  to 
adequately  reflect  existing  and  modified  (enlarged  chambers  and  conduits) 
conditions.  To  facilitate  the  meshing,  some  minor  simplifications  were 
incorporated  in  the  geometry  of  the  conduits  and  the  spillway  pier.  In  particular. 
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the  spillwaj'  pier  was  modeled  assuming  a  reduced  thickness  of  5  ft  (existing 
condition)  and  6  ft  (modified  condition),  instead  of  the  actual  thickness  of  8  ft. 
Both  the  mass  dens%  and  the  elastic  modulus  of  the  pier  concrete  were  adjusted 
by  a  fector  equal  to  the  ratio  betv^feen  the  actual  and  assumed  widths.  This 
approach  appropriately  captures  the  influence  of  the  pier  in  the  global  dynamic 
characteristics  of  the  section,  and  it  is  justified  in  those  cases  where  the  goal  of 
the  analysis  is  the  detemination  of  the  stresses  within  the  mass  concrete.  For  all 
cases  analyzed,  drj?  dommstream  conditions  were  assumed  (no  tailwater)  and  the 
reservoir  pool  elevation  was  assumed  at  gross  pool  elevation  (466.0  ft). 

Based  on  the  eigenvalue  analyses  of  the  3D  finite-element  models,  the  com¬ 
puted  fimdamental  vibration  period  of  monolith  14  (modified  condition)  on  rigid 
foundation  rock  with  empty  reservoir  was  0. 163  sec.  After  taking  into  account 
the  effects  of  water  and  foundation-rock  flexibility  according  to  the  simplified 
procedure,  the  resulting  value  of  the  fundamental  vibration  period  was  0.308  sec. 
For  monolith  21,  the  computed  fundamental  vibration  period  was  0. 183  sec, 
corresponding  to  rigid  foundation  rock  and  empty  reservoir  conditions.  This 
value  increased  to  0.261  sec  after  incorporating  water  and  foimdation-rock 
flexibilitj'  effects  according  to  the  simplified  procedure. 

Estimates  of  the  peak  seismic  responses  were  obtained  by  performing  static 
analysis  of  the  3D  models  subjected  to  the  tw'o  sets  of  earthquake-indueed  lateral 
forces  and  combining  the  corresponding  responses  using  the  square-root-of-the- 
sum-of-squares  (SRSS)  rule.  These  peak  values  were  then  combined  with  the 
static  effects  (dead  weight,  hydrostatic  pressure,  and  uplift)  to  provide  total 
response  estimates.  The  results  indicated  a  few  localized  areas  with  important 
tensile  demands.  For  example,  the  results  for  monolith  14  (modified  condition) 
showed  peak  vertical  tensile  stresses  within  the  apparent  dynamic  tensile  strength 
of  the  mass  concrete  (700  psi),  except  for  a  stress  concentration  of  1,143  psi  at 
the  upstream  heel  of  the  monolith  (corresponding  to  the  case  of  static  loads  plus 
horizontal  and  vertical  components  of  the  ground  motion).  This  stress  drops 
sharply  to  values  below  the  tensile  strength  within  a  distance  of  10  ft  up  above 
the  base.  The  results  for  monolith  21  also  indicated  peak  tensile  stresses  within 
the  apparent  djmamic  tensile  strength  of  the  mass  concrete,  excepting  some  areas 
of  stress  concentration.  For  example,  the  peak  vertical  stress  at  the  upstream  heel 
reached  893  psi  (corresponding  to  the  case  of  static  loads  plus  horizontal  and 
vertical  components  of  the  ground  motion),  but  the  value  quickly  drops  below  the 
tensile  strength  within  a  distance  of  10  ft  up  from  the  base.  In  spite  of  the  fact 
that  these  stress  concentrations  may  be  accentuated  by  the  assumed  fixed 
boundary  condition  at  the  dam  base,  high  tensile  stresses  are  very  likely  to  occur 
at  these  localized  regions.  However,  even  if  cracking  occurs  in  these  areas,  the 
extent  of  damage  is  expected  to  be  limited  and  this  is  still  considered  as 
acceptable  MCE  perfomance. 

Global  stability  anatyses  were  performed  for  monoliths  12  through  20.  The 
studies  included  t\^fo-dimensional  (2D)  conventional  stability?  analyses 
(determination  of  factor  of  safety  against  sliding  and  percentage  of  base  in 
compression),  conducted  in  compliance  with  EM  1 1 10-2-2200  (Headquarters, 
U.S.  Army  Corps  of  Engineers  (HQUSACE)  1995),  and  computation  of 
permanent  displacement  estimates,  following  the  procedure  indicated  in 
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EM  1110-2-6050  (HQUSACE  1999).  Two  possible  failure  surfaces  were 
considered:  a  concrete  lift  joint  near  the  base  of  the  dam  (lift  joint  failure),  and  a 
horizontal  plane  in  the  rock  mass  at  (or  just  below)  the  dam-fbundation  contact 
(near  contact  failure).  The  possibility  of  a  more  complex  failure  mode  occurring 
within  the  rock  mass  is  currently  under  investigation.  For  the  seismic  load  cases, 
the  reservoir  pool  elevation  was  assumed  at  446  ft.  This  value  was  based  on 
historical  data  and  it  represents  the  average  for  the  high  three  months  (May,  June, 
July).  Tailwater  level  was  assumed  at  elevation  160  ft.  Uplift  pressures  were 
computed  assuming  a  drain  efficiency  of  50  percent  for  all  monoliths 
investigated,  excepting  monolith  15,  for  which  60-percent  drain  efficiency  was 
assumed  based  on  actual  piezometer  data.  Each  monolith  was  analyzed 
considering  two  sets  of  shear  strength  parameters,  depending  on  the  assumed 
failure  surface.  For  lift  joint  failure,  a  cohesion  value  of  300  psi  and  a  friction 
angle  of  50  deg  were  selected  for  all  monoliths.  For  near  contact  failure,  zero 
cohesion  was  considered  and  the  friction  angle  was  computed  as  a  combination 
of  a  basic  rock-on-rock  friction  angle  (30  deg)  and  a  dilatancy  angle  representing 
localized  conditions  below  each  monolith.  Based  on  this  approach,  a  friction 
angle  of  55  deg  was  determined  for  the  analysis  of  near  contact  failure  in  all 
monoliths,  with  the  exception  of  monoliths  15  and  16,  which  were  assigned 
friction  angles  of  44  deg  and  53  deg,  respectively.  For  seismic  load  cases,  accept¬ 
able  values  of  factors  of  safety  were  foxmd  for  almost  all  the  monoliths  investi¬ 
gated,  with  the  exception  of  monolith  15.  The  computed  sliding  factor  of  safety 
under  MCE  conditions  was  1.0,  and  therefore  this  monolith  did  not  meet  the  cor¬ 
responding  MCE  sliding  stability  criteria.  However,  based  on  the  computed  esti¬ 
mates  of  sliding  permanent  displacement  (less  than  1  in.),  its  performance  was 
considered  acceptable. 
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2  Selection  of  Ground 
Motions 


Maximum  Credible  Earthquake 

The  earthquake  excitation  was  specified  by  the  Sacramento  District 
according  to  the  most  recent  seismic  hazard  studies.  These  studies  defined  the 
controlling  MCE  as  an  event  of  magnitude  6.5  at  a  source-to-site  distance  of 
14  km,  on  the  eastern  branch  of  the  Bear  Mountains  fault  zone.  Based  on  this,  the 
characterization  of  the  corresponding  site  ground  motions  was  performed  by 
considering  four  different  attenuation  relationships.  The  peak  horizontal  ground 
acceleration  values  corresponding  to  the  median  (50'*'  percentile)  and  median  plus 
one  standard  deviation  (S^***  percentile)  were  determined  as  0.24  g  and  0.38  g, 
respectively  (URS  Corporation  2001).  Based  on  the  magnitude  and  distance  of 
the  MCE,  it  was  concluded  that  faults  directivity  effects  were  not  significant  and 
the  same  response  spectram  was  specified  for  the  two  horizontal  components. 

The  corresponding  vertical  response  spectrum  was  initially  defined  by  uni¬ 
form  scaling  of  the  MCE  horizontal  spectrum  using  a  reduction  factor  of  2/3 
(URS  Corporation  2002).  Therefore,  the  vertical  response  spectrum  was  deter¬ 
mined  based  on  a  constant  factor  for  all  periods.  This  recommendation  was  later 
modified  and  the  resulting  vertical  response  spectram  was  defined  using  a 
period-dependent  scaling  factor  (URS  Corporation  2003).  The  resulting 
5-percent-damping  response  spectra  for  horizontal  and  vertical  motion  are  showTi 
in  Figure  1. 


Ground  Motion  Time  Histories 

Three  sets  of  earthquake  time  histories,  each  consisting  of  three  components, 
w'ere  developed  for  the  MCE  (URS  Corporation  2002, 2003).  These  time  histories 
were  generated  by  modifying  a  suite  of  recorded  ground  motioiK  so  their  response 
spectra  would  match  the  corresponding  5 -percent-damping  response  spectra 
showTi  in  Figure  1.  The  earthquAe  records  used  to  develop  the  spectrally  matched 
time  histories  are  presented  in  Table  1 .  Out  of  the  two  horizontal  components  for 
each  time  history',  tire  one  associated  with  the  largest  peak  ground  velocity  (PGV) 
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Period  [sec] 


Figure  1.  MCE  5-percent-damped  horizontal  and  vertical  response  spectra 


Table  1 

Input  Ground  Motions  Used  in  the  Analyses 

Earthquake 

Recorded  Ground  Motions 

Modified  Time  Histories 

Mw 

Station 

KSSI 

■1 

IH 

Direction 

1971  San  Fernando 

6.6 

Pasadena  -  Old 
Seism.  Lab. 

19 

180 

0.09 

0,38 

27.0 

Cross  Ch. 

270 

0.20 

0.38 

34.8 

Us/Ds 

Vertical 

0.09 

0.30 

13.5 

Vertical 

1979  Imperial  Valley 

6.5 

Cerro  Prieto 

26 

147 

0.17 

0.38 

23,8 

Us/Ds 

237 

0.16 

0.38 

23.1 

Cross  Ch. 

Vertical 

0.21 

0.33 

11.5 

Vertical 

1986  Chalfant  Valley 

6.2 

Bishop  -  Paradise 
Lodge 

23 

70 

0.16 

0.38 

28.8 

Cross  Ch. 

160 

0,16 

0.38 

29.4 

Us/Ds 

Vertical 

0,13 

0.31 

11.7 

Vertical 

was  selected  as  the  horizontal  upstream-downstream  component  to  be  used  in  the 
two-dimensional  analyses.  Figures  2  through  4  show  the  selected  horizontal  and 
vertical  components  for  the  spectrally  matched  time  histories. 

The  initial  phase  in  the  evaluation  of  the  seismic  performance  of  a  concrete 
gravity  section  using  linear  time-history  analysis  is  primarily  concerned  with  the 
determination  of  the  peak  values  of  relevant  response  quantities,  such  as  the  max¬ 
imum  principal  tensile  stresses,  and  their  spatial  distribution.  For  the  purpose  of 
determination  of  peak  responses,  the  same  duration  of  analysis  (30  sec)  was  used 
for  the  three  sets  of  input  motions  considered  in  this  study.  The  decision  to  use  a 
uniform  duration  for  all  cases  was  made  based  on  preliminary  analyses  that 
showed  that  this  duration  was  adequate  to  bracket  the  peak  values  and  capture  the 
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Figure  4.  1986  Chalfant  Valley  earthquake  spectrally  matched  time  histories 


relevant  features  of  the  response.  Table  2  shows  a  comparison  of  the  duration 
parameters  corresponding  to  the  horizontal  components  of  the  selected  input 
ground  motions  for  the  original  records  and  for  the  30-sec  truncated  records.  Fig¬ 
ures  5  and  6  highlight  the  close  match  between  the  individual  response  spectra 
for  the  30-sec  truncated  input  time  histories  and  the  corresponding  targets  for  the 
horizontal  and  vertical  directions,  respectively. 


Table  2 

Duration  Parameters  for  Horizontal  Components  of  Input  Ground  Motions 

Parameter 

Original  records 

Truncated  records  (30  fseci)  1 

San  Fernando 

Imperial 

Valley 

Chalfant 

San  Fernando 

Imperial 

Valley 

Chalfant 

Total  Arias  Intensity  4,  in/sec 

103,15 

321.02 

140.17 

103.15 

265.33 

138.16 

Starting  Time  (5  percent  4),  sec 

2.24 

4.88 

1.74 

2.24 

4.54 

1.71 

Ending  Time  (95  percent  4),  sec 

16.30 

44.36 

17.35 

16.30 

27.40 

16.08 

Arias  Intensity  Duration  Tiafs-gs),  sec 

14.06 

39.48 

15.61 

14.06 

22.86 

14.36 

Bracketed  Duration  T0.05,  sec 

22.35 

57.37 

35.12 

22.35 

28.53 

29.53 
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re  5.  Comparison  of  5-percent-damped  horizontal  response  spectra  for  truncated  time 
histories 


Figure  6,  Comparison  of  5-percent-damped  vertical  response  spectra  for  truncated  time  histories 
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3  Dynamic  Stress  Analyses 


Introduction 

This  section  presents  the  results  of  the  dynamic  stress  analyses  performed  on 
two-dimensional  (2D)  finite-element  models  of  monoliths  14  and  21,  subject  to 
ground  motion  time  histories  that  are  representative  of  the  MCE.  Different  analy¬ 
ses  are  performed  using  the  computer  program  EAGD-84.  The  objective  of  these 
analyses  is  not  only  to  quantify  the  seismic  performance  of  the  monoliths  by 
determining  the  stress  responses,  but  also  to  investigate  the  sensitivity  of  the 
dynamic  response  characteristics  and  corresponding  results  with  respect  to  vari¬ 
ations  in  the  modeling  assumptions. 


Material  Properties 

A  number  of  field  and  laboratory  investigations  have  been  performed  over 
the  years  to  characterize  the  material  parameters  of  the  mass  concrete  and  the 
foundation  rock.  Table  3  suimnarizes  some  of  the  relevant  material  properties  for 
the  mass  concrete,  as  defined  by  the  Sacramento  District.  Because  it  is  difficult  to 
accurately  define  the  thickness  of  the  layer  of  rich  concrete,  it  is  assumed  that  the 
mass  concrete  is  characterized  by  the  material  properties  corresponding  to  the 
lean  concrete  mix. 


Table  3 

Material  Parameters  for  Mass  Concrete 

Material  Property 

Value 

Unit  weight  ,  pcf 

158 

Concrete  compressive  strength  f^' 

4,375 

Splitting  tensile  strength  ,  psi 

360 

^  st 

Static  modulus  of  rupture  ,  psi 

540 

Dynamic  modulus  of  rupture  ,  psi 

810 

Static  modulus  of  elasticity  ,  psi 

3.6. 10® 

Dynamic  modulus  of  elasticity  ,  psi 

5.9. 10® 

Poisson’s  ratio 

0.19 
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The  values  indicated  above  for  the  unit  weight,  modulus  of  elasticity,  and 
Poisson’s  ratio  are  the  same  as  those  used  in  a  previous  seismic  evaluation  study 
(Hall  et  al.  1989)  as  well  as  in  the  stress  analyses  conducted  by  the  Sacramento 
District  using  Chopra’s  simplified  procedure. 

The  tensile  behavior  of  the  mass  concrete  plays  a  critical  role  in  the 
evaluation  of  the  seismic  performance  of  a  gravity  section,  and  it  is  therefore 
important  to  understand  the  basis  for  the  strength  value  indicated  in  Table  3.  The 
1989  study  adopted  a  value  of  700  psi  for  the  apparent  tensile  strength  of  the  lean 
concrete  (this  parameter  is  conceptually  equivalent  to  the  dynamic  modulus  of 
rapture).  This  was  based  on  the  analysis  of  a  series  of  static  and  djmamic  split 
cylinder  tests  conducted  by  the  U.S.  Bureau  of  Reclamation  in  1983  and  the 
University  of  California  at  Berkeley  in  1986.  In  2001,  the  U.S.  Bureau  of 
Reclamation  performed  a  new  series  of  split  cylinder  tests  on  12-in.-diam 
concrete  specimens.  The  corresponding  results  indicated  an  average  splitting 
tensile  strength  (^)  of  360  psi  for  lean  concrete  (with  test  values  that  ranged 
from  216  to  428  psi).  The  value  of  the  static  modulus  of  rapture  was  computed 
as  //*  =  1 .5^ ,  and  a  dynamic  increase  of  50  percent  was  assumed  for  the 
djTiamic  modulus  of  rupture,  that  is/,.  =1.5  //  =  2.25  / .  This  leads  to  the  value 
of  8 10  psi  indicated  in  Table  2.  However,'  a  lower  value  of  tensile  strength 
(700  psi)  was  considered  as  the  maximum  allowable  tensile  stress  for  evaluation 
of  the  dam  performance  in  the  stress  analyses  conducted  by  the  Sacramento 
District. 

Table  4  summarizes  the  relevant  material  properties  for  the  foundation  rock, 
as  defined  by  the  Sacramento  District.  These  values  correspond  to  the  average 
values  indicated  in  the  1989  study  (Hall  et  al.  1989),  and  they  were  based  on  the 
results  of  a  series  of  geophysical  surveys  conducted  by  the  Corps  of  Engineers  in 
1989.  These  investigations  focused  mainly  on  the  left  and  right  abutments  of  the 
dam,  and  it  was  recognized  that,  in  some  cases,  the  geophysical  studies  did  not 
penetrate  deep  enough  to  substantially  engage  the  foundation  rock  in  all  locations 
tested.  In  particular,  the  values  of  djmamic  modulus  of  elasticity  ranged  from 
5.8  •  10®  psi  to  1 1.0  •  10®  psi,  and  the  value  of  7.9  •  10®  psi  has  been  accepted  by 
the  Sacramento  District  as  a  reasonable  estimate  for  the  dynamic  modulus  of 
elasticity  of  the  foundation  rock  at  Folsom  Dam. 


Table  4 

Material  Parameters  for  Foundation  Rock 

Material  Property 

Value 

Unit  weight  Jj- ,  pcf 

171 

Dynamic  modulus  of  ela^icity  ,  psi 

7.9. 10® 

Poisson's  ratio  Vj- 

0.25 
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EAGD-84  Analyses 

Finite-element  models 

The  EAGD-84  computer  program  was  developed  at  the  University  of  Cali¬ 
fornia  at  Berkeley  to  evaluate  the  seismic  response  of  two-dimensional  sections 
of  concrete  gravity  dams,  including  the  effects  of  dam-water  interaction,  dam- 
foundation  rock  interaction,  and  energy  absorption  effects  at  the  bottom  of  the 
reservoir  (Fenves  and  Chopra  1984). 

The  general  analytical  model  is  based  on  the  substructure  method,  and  it 
assumes  linear  behavior  for  the  dam-water-foundation  system.  The  resulting 
equations  of  motion  are  solved  in  the  frequency  domain.  The  model  assumes  a 
concrete  section  with  a  horizontal  dam-foundation  interface,  and  constant  hyste- 
retic  damping  is  used  to  represent  energy  dissipation  in  the  mass  concrete.  The 
foundation  rock  underlying  the  dam  is  idealized  as  a  homogeneous,  isotropic, 
visco-elastic  half-plane.  The  water  impounded  in  the  reservoir  is  modeled  as  a 
fluid  domain  of  constant  depth  and  infinite  length  along  the  upstream  direction, 
and  the  energy  absorption  associated  with  reservoir  bottom  materials  is 
quantified  by  a  wave  reflection  coefficient.  The  earthquake  excitation  is  defined 
by  two  components  (upstream/downstream  and  vertical)  of  the  ffee-field  ground 
acceleration,  and  the  corresponding  ground  acceleration  time-histories  are 
applied  at  the  base  of  the  dam. 

Geometry  and  relevant  dimensions  of  monoliths  21  (nonoverflow)  and  14 
(overflow)  are  presented  in  Figures  7  and  8.  The  middle  cross  sections  of  these 
monoliths  were  idealized  by  two-dimensional  finite-element  models  consisting  of 
four-node  quadrilateral  elements.  The  same  basic  set  of  material  properties  was 
uniformly  assigned  to  the  mass  concrete,  and  therefore  monolith  2 1  was  modeled 
with  only  one  material  assigned  to  the  whole  section.  The  cross  section  of  this 
monolith  was  discretized  with  1712  quadrilateral  elements,  as  shown  in  Figure  9. 
It  is  important  to  notice  that  the  mesh  exhibits  only  eight  elements  along  the  dam- 
foundation  interface.  As  a  result  of  existing  limitations  in  the  implementation  of 
the  EAGD-84  model,  this  transition  to  a  reduced  number  of  elements  along  the 
dam-foundation  contact  is  necessary  to  guarantee  that  the  foundation  rock  is 
represented  by  a  dynamic  stiffiiess  matrix  that  is  fully  consistent  with  the 
assumption  of  a  visco-elastic  half  plane.  The  reduced  number  of  elements  at  the 
base  of  the  model  will  require  a  careful  interpretation  of  the  corresponding  stress 
results.  In  spite  of  the  fact  that  this  number  of  elements  should  be  sufficient  to 
capture  the  general  characteristics  of  the  stress  distribution  at  the  base  of  the 
stmcture,  this  modeling  limitation  will  produce  underestimation  of  the  stresses  in 
those  areas  with  high  stress  gradients  (such  as  the  heels).  Therefore,  it  will  be 
important  to  evaluate  the  stress  results  not  only  in  terms  of  their  peak  values  but 
also  in  terms  of  their  distribution.  This  must  be  also  considered  when  conducting 
the  comparison  of  these  results  with  those  obtained  by  the  Sacramento  District 
using  refined  3D  models  of  the  monoliths. 
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Figure  7.  Geometry  of  nonoverflow  monolith  21 
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Figure  9.  Finite-element  mesh  for  nonoverflow 
monolith  21 

The  middle  cross  section  of  monolith  14  was  discretized  with  1936  quadri¬ 
lateral  elements,  as  shown  in  Figures  10  through  12.  As  in  the  previous  case,  it 
becomes  necessary  to  incorporate  a  transition  to  a  reduced  number  of  elements 
along  the  base.  In  addition,  since  this  is  an  overflow  section,  the  presence  of  the 
spillway  pier  and  the  openings  corresponding  to  the  upper  and  lower  chambers 
and  conduits  must  be  accounted  for  in  the  model.  The  two-dimensional  finite- 
element  models  used  for  these  studies  are  based  on  the  assumption  that  the 
section  is  regular  along  the  direction  normal  to  its  plane,  i.e.,  no  changes  in 
material  properties  or  dimensions.  Therefore,  in  this  case  it  is  necessary  to 
account  for  the  local  variations  in  thickness  of  the  section  by  appropriately 
modifying  the  corresponding  material  properties.  The  motivation  for  this 
approach  is  to  represent-globally-the  effects  that  these  irregularities  have  on  the 
dynamic  characteristics  of  the  section,  and  it  is  achieved  by  introducing 
equivalent  mass  and  stiffiiess  distributions.  Of  course,  a  comprehensive 
investigation  of  the  stress  concentration  effects  that  are  induced  by  the  presence 
of  these  thickness  irregularities  would  require  three-dimensional  models  such  as 
those  used  by  the  Sacramento  District  in  previous  studies.  Three  different  finite- 
element  models  were  developed  for  monolith  14:  (a)  the  first  one  accounting  only 
for  the  presence  of  the  spillway  pier  (thus  ignoring  the  voids  created  by  the 
chambers  and  conduits),  (b)  the  second  one  considering  not  only  the  spillway 
pier  but  also  the  chambers  and  conduits  by  means  of  equivalent  material 
properties,  and  (c)  the  third  model  incorporates  the  voids  corresponding  to  the 
chambers  while  still  using  equivalent  material  properties  for  the  pier  and  the 
conduits.  Table  5  summarizes  the  material  properties  used  in  the  finite-element 
models  for  monoliths  14  and  21. 
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Figure  1 0.  Finite-element  mesh  for  spillway  monolith 

14  (Model!) 


Figure  1 1 .  Finite-element  mesh  for  spillway  monolith  1 4 
(Model  2) 
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Figure  12.  Finite-element  mesh  for  spiliway  monolith 
14  (Model  3) 


Table  5 

Material  Parameters  Used  in  Finite-Element  Models 


Monoliths  14 

Monolith  14  || 

Material  Property 

and  21  (Mass 
Concrete) 

Zone  A 
(Pier) 

Zone  B 
(Chambers) 

Zone  C 
(Conduits) 

Unit  weight  Yf  .  pcf 

158 

25.3 

102.7 

119.1 

Dynamic  modulus  of  elasticity  Ef ,  psi 

5.9-10® 

0.94-10® 

3.84-10® 

4.45-10® 

Poisson’s  ratio  Vf 

0.19 

0.19 

0.19 

0.19 

The  dynamic  stress  analyses  described  in  the  next  sections  included  the 
effects  of  static  loads  (self  weight  and  external  hydrostatic  loads),  and  the  reser¬ 
voir  water  surface  was  assumed  at  gross  pool  elevation  (466  ft).  However,  effects 
of  uplift  pressures  or  tailwater  were  not  included  in  the  stress  computations  pre¬ 
sented  in  this  report. 

Analysis  parameters 

In  addition  to  the  standard  set  of  material  properties  described  in  the  previous 
section,  the  formulation  implemented  in  the  EAGD-84  computer  program 
requires  the  definition  of  a  series  of  additional  parameters  that  are  not  typically 
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required  when  using  multi-purpose  analysis  programs  such  as  SAP2000.  The 
most  relevant  of  these  parameters  is  the  wave  reflection  coefficient  (  a  ).  This 
coefficient  represente  the  dissipation  of  hydrodynamic  pressure  waves  in  the 
reservoir  bottom.  When  pressure  waves  strike  the  bottom  surface,  part  of  the 
propagating  energy  is  reflected  back  into  the  reservoir  and  another  part  is 
transmitted  to  die  foundation  medium  because  of  its  finite  stiffiiess.  This  eneigj' 
may  be  assumed  to  continue  propagating  towards  infinity,  actually  constituting 
an  equivalent  energy-dissipating  mechanism  introducing  radiation  damping  and 
effectively  reducing  die  system’s  response.  The  wave  reflection  coefficient  is 
defined  as  follows: 


l-qC 

a  = - - — 

\  +  qC 

where  C  is  velocity  of  pressure  waves  in  water,  and  q  is  given  by 


(1) 


In  this  equation,  p,^is  the  mass  density  of  water,  tiiat  is,  where 

indicates  the  unit  weight  and  g  is  the  acceleration  of  gravity.  The  parameters 
p^  and  denote  the  mass  density  and  elastic  modulus  of  the  reservoir  bottom 
material,  respectively.  Assuming  values  of  4,730  ft/sec  and  62.4  Ib/ft^  for  the 
velocity  of  pressure  waves  and  imit  weight  of  the  water,  and  assuming  that  the 
properties  of  the  reservoir  bottom  material  are  the  same  as  those  adopted  for  the 
foundation  rock  {E^=Ef,y^=jf)as  listed  in  Table  3,  a  value  of  0.79  is 

obtained  for  the  wave  reflection  coefficient.  This  value  is  representative  of 
energy  absorption  in  the  rock  material  itself,  and  it  does  not  actually  account  for 
any  additional  absorption  due  to  actual  sediment  materials  deposited  at  the 
bottom  of  the  reservoir.  It  is  important  to  highlight  that  results  from  field 
measurements  conducted  in  1995  indicated  the  presence  of  sediment  layers  with 
a  significant  thickness. 

A  value  of  0.75  was  adopted  for  the  wave  reflection  coefficient  in  previous 
seismic  evaluation  studies  conducted  by  the  Sacramento  District.  This  value  is 
also  adopted  for  the  stress  analyses  presented  in  this  report.  Figure  13  shows  the 
variation  of  the  wave  reflection  coefficient  as  a  function  of  the  elastic  modulus  of 
the  reservoir  bottom  material,  E^ ,  which  is  represented  in  terms  of  the 
dimensionless  ratio  E^l  Ej- .  The  wave  reflection  coefficient  is  presented  for 
different  values  of  the  mass  density'  of  the  reservoir  bottom  materials,  which  are 
expressed  in  terms  of  the  rmit  weight  ratio  y^ty^.  Based  on  the  figure,  it  can  be 
concluded  that  the  value  adopted  for  the  previous  seismic  evaluation  studies 
reasonably  represents  a  range  of  possible  combinations  of  actual  values  of  E^ 
and  p,.  at  the  site. 
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Ratio  E/E, 


Figure  13.  Variation  of  the  wave  reflection  coefficient  as  a  function  of  the  material  parameters  of 
the  reservoir  bottom  material 

In  addition  to  the  wave  reflection  coefficient,  it  is  required  to  quantify  two 
other  parameters  representing  different  sources  of  damping  considered  in  the 
analysis.  One  is  a  constant  hysteretic  damping  model  for  the  dam  concrete  {r]^) 
and  the  other  is  a  constant  hysteretic  damping  model  for  the  foundation  rock 
medium  ( /^^ ).  A  value  of  hysteretic  damping  coefficient  of  10  percent 

(77^  =0.10)  is  typically  recommended  for  the  dam,  and  this  value  corresponds  to 
a  5  percent  viscous  damping  ratio  in  all  vibration  modes  of  the  dam.  This  value  is 
considered  suitable  for  the  large  motions  and  associated  high  stresses  expected 
during  strong  earthquake  shaking.  A  similar  level  of  hysteretic  damping  is 
considered  for  the  foundation  rock  ( 7?^^  =  0. 10 ),  following  the  same  assumptions 

employed  for  the  previous  studies. 


Results  -  monolith  21 

Tables  6  and  7  show  the  first  ten  natural  frequencies  for  monolith  21,  corre¬ 
sponding  to  empty  reservoir  conditions  and  two  different  foundation  conditions. 
The  natural  frequencies  were  calculated  for  rigid  {Ef=<x>)  and  flexible 

{Ef -1.9  •  10®  psi)  foundation  cases.  The  incorporation  of  foundation  flexibility 

causes  an  increase  of  about  1 8  percent  in  the  fundamental  period  of  the  dam- 
foundation  system  with  respect  to  the  dam  on  rigid  base.  The  influence  of  the 
flexibility  increases  all  the  natural  periods,  with  a  significant  effect  on  the  period 
corresponding  to  the  third  vibration  mode.  This  mode  is  associated  with  a 
predominantly  vertical  response  of  the  monolith. 
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Table  6 

Natural  Frequencies  -  Monolith  21  with  Rigid  Foundation  (Empty 
Reservoir) 

— II II 1 1  |i|  1 1  II 1 

Frequency 

Period 

[Hz] 

[sec] 

1 

5,44 

0.184 

2 

12.09 

0.083 

3 

16.93 

0.059 

4 

22.53 

0.044 

5 

34.20 

0.029 

6 

36.30 

0.028 

7 

47.54 

0.021 

8 

50.11 

0.020 

9 

52.33 

0.019 

10 

57.77 

0.017 

Table  7 

Natural  Frequencies  -  Monolith  21  with  Flexible  Foundation 
(Empty  Reservoir) 

Mode 

Flexible  Foundation  I 

Frequency 

Period 

[Hz] 

[sec] 

1 

4.62 

0.217 

2 

9.12 

0,110 

3 

9.19 

0.109 

4 

16.52 

0.061 

5 

26.66 

0.038 

6 

27.59 

0.036 

7 

39.83 

0.025 

8 

42.74 

0.023 

9 

45.01 

0.022 

10 

51 .24 

0.020 

Figures  14  and  15  show  the  location  of  the  first  five  natural  periods  with 
respect  to  the  MCE  5 -percent-damped  horizontal  and  vertical  response  spectra. 
Figure  13  corresponds  to  rigid  foundation  conditions,  whereas  Figure  14  repre¬ 
sents  the  effects  of  foundation  flexibility.  It  is  important  to  mention  that  these 
natural  periods  correspond  to  empty  reservoir  conditions  (as  computed  by 
EAGD-84),  and  that  the  actual  dynamic  response  characteristics  of  the  system 
will  be  affected  by  the  corresponding  hydrodynamic  effects. 

Figure  16  show's  the  distribution  of  static  normal  vertical  stress  (Sy)  along  the 
dam-foundation  interface,  considering  only  self-weight  and  hj'drostatic  loading. 
The  results  are  presented  for  both  rigid  and  flexible  foundation  conditions.  The 
standard  stress  output  in  EAGD-84  corresponds  to  the  element  centers,  and  these 
values  w'ere  averaged  between  adjacent  elements  along  the  dam-foundation  inter¬ 
face  to  compute  the  corresponding  nodal  values.  The  stress  values  for  tiie  heel  and 
toe  nodes  were  computed  by  linear  extrapolation.  The  figure  shows  tiiat  for  the 
assumed  reservoir  level  (gross  pool  elevation  at  466  ft),  the  static  normal  stress  at 
the  heel  of  the  dam  teaches  values  of  60  and  -65  psi,  for  rigid  and  flexible  base 
conditions,  respectively.  The  distribution  of  the  normal  stress  Sy  across  the  section 
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Acceleration  [g]  ^  Acceleration  [g] 


!  14.  Natural  periods  and  MCE  spectra  -  monolith  21  with  rigid  foundation 


Figure  15.  Natural  periods  and  MCE  spectra  -  monolith  21  with  flexible  foundation 
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is  depicted  by  the  contour  plots  shown  in  Figures  17  and  18,  for  rigid  and  flexible 
foiuidation  conditions,  respectively. 

Since  the  structure  is  not  symmetric,  the  dynamic  stresses  along  the  upstream 
and  downstream  faces  depend  on  the  sign  convention  used  to  scale  the  ground 
acceleration  records.  To  identify  the  critical  loading  scenario,  several  analyses 
were  performed  for  each  earthquake,  combining  the  initial  static  loading  with 
each  of  the  two  possible  directions  of  the  horizontal  ([+H],  [-H])  and  vertical 
([+V],  [-V])  components.  This  resulted  in  four  cases  for  each  earthquake: 

[+H+V],  [+H-V],  [-H+V],  and  [-H-V].  A  total  of  12  different  cases  were  obtained 
when  these  combinations  of  horizontal  and  vertical  components  were  considered 
for  each  one  of  the  three  earthquakes.  Furthermore,  to  investigate  the  influence  of 
the  vertical  component  of  the  seismic  input,  dynamic  stress  analyses  were 
performed  considering  only  the  horizontal  component  of  each  ground  motion. 
Considering  the  two  possible  directions  for  the  horizontal  component  of  each 
earthquake,  this  resulted  in  six  additional  ground  motion  cases.  Since  the  damage 
to  a  concrete  gravity  section  is  essentially  associated  with  the  development  of 
tensile  stresses,  the  critical  loading  scenario  for  each  earthquake  was  identified  as 
the  combination  that  generated  the  highest  tensile  stress  value  in  the  section. 

Tables  8, 9,  and  10  summarize  some  of  the  results  from  the  dynamic  stress 
analyses  of  monolith  21 .  For  each  earthquake,  the  tables  show  the  peak  values  of 
maximum  principal  stresses,  their  location,  and  time  of  occurrence.  The  stress 
values  represent  total  stresses  as  they  include  the  static  loading  effects.  The  pur¬ 
pose  of  Tables  8  through  10  is  to  identify  the  critical  component  combination. 

The  reported  principal  stress  values  are  element  centroidal  stresses  and  therefore 
their  location  corresponds  to  the  centroid  of  the  corresponding  element. 
Additional  post-processing  was  performed  in  the  output  data  to  generate  the 
nodal  stress  values  used  in  the  contour  plots  shown  later. 
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1 7.  Static  stress  Sy  -  monolith  21  with  rigid 
foundation 


1 8.  Static  stress  Sy  -  monolith  21  with  flexible 
foundation 


Peak  Values  of  Maximum  Principal  Stresses  -  Monolith  21 
San  Fernando  Earthquake 


on 


ase  (Heel) 


Upstream 


Downstream 


Base  (Heel) 


Upstream 


Downstream 


ase  (Heel) 


Upstream 


Downstream 


Base  (Heel) 


Upstream 


Downstream 


Base  (Heel) 


Upstream 


Downstream 


Base  (Heel) 


Upstream 


Downstream 


4.85 

8.75 

20.53 

196.31 

61.87 

196.31 

4.85 

8.75 

20.53 

196.31 

63.64 

192.92 

4.85 

8.75 

20.53 

196.31 

61.87 

196.31 

4.85 

8.75 

20.53 

196.31 

63.64 

192.92 

4.85 

8.75 

20.53 

196.31 

61.87 

196.31 

4.85 

8.75 

20.53 

196.31 

60.45 

199.25 

Table  9 

Peak  Values  of  Maximum  Principal  Stresses  -  Monolith  21 
Imperial  Valley  Earthquake 


Location 


Base  (Heel) 


Upstream 


Downstream 


Base  (Heel) 


Upstream 


Downstream 


Base  (Heel) 


Upstream 


Downstream 


Base  (Heel) 


Upstream 


Downstream 


ase  (Heel) 


Upstream 


ownstream 


ase  (Heel) 


earn 


earn 


4.85 

8,75 

19.6 

633 

20.53 

196,31 

21.7 

477 

61.87 

196.31 

21.8 

504 

4.85 

8.75 

24,7 

568 

20.53 

196.31 

21.8 

513 

63.64 

192.92 

21.7 

492 

4.85 

8.75 

19.6 

704 

20.53 

196.31 

22.2 

674 

63.64 

192.92 

12.2 

723 

4.85 

8,75 

21.7 

661 

20.53 

196,31 

15.3 

598 

63.64 

192.92 

15.2 

524 

4.85 

8.75 

5.9 

616 

20.53 

196.31 

21.8 

499 

61.87 

196.31 

15.3 

575 

4.85 

8.75 

22.1 

737 

20.53 

12,2 

705 

61.87 

196.31 

22.2 

691 
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Table  10 

Peak  Values  of  Maximum  Principal  Stresses  -  Monolith  21  - 
Chalfant  Earthquake 

Case 

Location 

X 

[ft] 

Y 

[ft] 

Time 

[sec] 

CTmax  Foundation 
[psi] 

+H 

Base  (Heel) 

4.85 

8.75 

6.0 

699 

+H 

Upstream 

20.53 

196.31 

7.0 

558 

+H 

Downstream 

61.87 

196.31 

2.7 

645 

-H 

Base  (Heel) 

4.85 

8.75 

2.7 

629 

-H 

Upstream 

20.53 

196.31 

2.7 

635 

-H 

Downstream 

61.87 

196.31 

7.0 

562 

+H+V 

Base  (Heel) 

4.85 

v;  8.75 

6.0 

749 

+H+V 

Upstream 

20.53 

196.31 

10.0 

558 

+H+V 

Downstream 

61.87 

196.31 

2.7 

719 

+H-V 

Base  (Heel) 

4.85 

8.75 

6.0 

665 

+H-V 

Upstream 

20.53 

196.31 

7.0 

682 

+H-V 

Downstream 

61.87 

196.31 

2.7 

584 

-H+V 

Base  (Heel) 

4.85 

8.75 

2.7 

602 

-H+V 

Upstream 

20.53 

196.31 

2.7 

596 

-H+V 

Downstream 

61.87 

196.31 

7.0 

678 

-H-V 

Base  (Heel) 

4.85 

8.75 

2.7 

656 

-H-V 

Upstream 

20.53 

196.31 

2.7 

678 

-H-V 

Downstream 

60.45 

199.25 

7.0 

563 

For  a  given  earthquake,  three  different  stress  values  are  reported  for  each 
combination  of  components,  and  they  are  the  peak  centroidal  stresses  along  the 
base,  along  the  upstream  face,  and  along  the  downstream  face.  For  the  San  Fer¬ 
nando  earthquake  (Table  8),  the  peak  stress  value  was  found  at  the  heel  of  the 
dam  and  it  reached  a  value  of  757  psi  for  the  case  [+H-V].  This  stress  value 
corresponds  to  the  center  of  the  comer  element  at  the  heel  of  the  dam,  located  at 
4.85  ft  from  the  upstream  face  and  8.75  ft  above  the  base.  Although  this  comer 
region  exhibits  the  highest  tensile  stresses  in  the  section,  the  corresponding  stress 
gradient  varies  very  rapidly  and  the  values  are  significantly  reduced  towards  the 
interior  of  the  section.  The  same  component  combination  generated  the  highest 
tensile  stress  (666  psi)  along  the  upstream  face,  at  a  location  196  ft  above  the 
base.  The  highest  tensile  stress  along  the  downstream  face  was  found  for  the  case 
[-H+V]  and  it  reached  674  psi,  at  a  location  196  ft  above  the  face  near  the 
change  of  slope.  The  results  for  the  Imperial  Valley  earthquake  (Table  9) 
indicated  a  lower  tensile  stress  at  the  heel  (737  psi),  but  higher  peak  stresses 
along  both  downstream  and  upstream  faces  (723  and  705  psi,  respectively)  when 
compared  to  the  previous  results  obtained  for  the  San  Fernando  earthquake.  The 
peak  stresses  at  the  heel  and  along  the  upstream  face  occurred  for  the  case  [-H- 
V],  whereas  the  peak  stress  along  the  downstream  face  corresponded  to  the  case 
[+H+V].  Another  significant  difference  with  respect  to  the  previous  case  is  given 
by  the  times  of  occurrence,  since  the  peak  values  occurred  much  later  into  the 
analyses.  The  results  for  the  Chalfant  earthquake  (Table  10)  revealed  apeak 
tensile  stress  at  the  heel  of  749  psi  for  the  case  [+H+V].  The  critical  tensile 
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stresses  along  tiie  downstream  and  upstream  faces  reached  719  psi  (case  [+H+V]) 
and  682  psi  (c^e  [+H-V]),  respectively,  and  they  occurred  at  the  same  elevations 
as  in  the  previous  cases. 

In  all  three  earthquakes,  the  critical  tensile  stress  conditions  were  reached  for 
cases  involving  the  combination  of  both  horizontal  and  vertical  ground  motion 
components.  The  reported  peak  stresses  and  their  times  of  occurrence  depend  on 
the  particular  earthquake  considered,  but  in  general,  the  peak  values  exhibit 
roughly  the  same  magnitude.  The  critical  combination  conditions  are  controlled 
in  all  cases  by  the  stress  at  the  heel  and  they  are  given  by  [+H-V],  [-H-V],  and 
[+H+V]  for  the  San  Fernando,  Imperial  Valley,  and  Chalfant  earthquakes, 
respectively. 

Figures  19  through  21  show  the  time  historj'  variation  of  the  maximum 
principal  (centroidal)  stresses  at  the  base,  downstream  face,  and  upstream  fece, 
respectivety,  for  the  critical  component  combination  corresponding  to  the  San 
Fernando  earthquake  (case  [+H-V]).  These  principal  stress  time  histories  reach 
peak  values  of  757,  641,  and  666  psi,  respectively.  In  general,  the  time  histories 
show'  that  the  most  significant  tensile  demands  clearly  occur  at  isolated  instances. 
Figures  22  through  24  show'  the  time  historj'  variation  of  the  maximum  principal 
stresses  at  the  base,  dow'nstream  face,  and  upstream  face,  respectively,  for  the 
critical  component  combination  corresponding  to  the  Imperial  Valley  earthquake 
(case  [-H-V]).  These  principal  stress  time  histories  reach  peak  values  of  737, 691, 
and  705  psi,  respectively.  Figures  25  through  27  show  similar  time  histories  for 
the  Chalfant  earthquake  (case  [+H+V).  The  stress  histories  reach  peak  values  of 
745,  719,  and  558  psi,  respectively.  As  mentioned  before,  the  maximum  principal 
stress  at  the  base  (heel)  show's  the  largest  tensile  peak  within  the  section  for  the 
three  ground  motions  considered. 

The  distribution  of  the  maximum  tensile  demands  across  the  section  is  typi¬ 
cally  depicted  by  means  of  contour  plots  that  indicate  the  peak  values  reached  by 
a  given  stress  quantity  over  the  duration  of  the  analysis.  Figures  28  through  30 
show'  the  contour  plots  of  the  maximum  values  of  normal  vertical  stresses  (Sj.)  for 
the  critical  component  combination  corresponding  to  the  San  Fernando,  Imperial 
Valley,  and  Chalfant  earthquakes,  respectively.  As  in  the  static  case,  normal 
vertical  stresses  were  selected  for  these  contour  plots  to  facilitate  the  comparison 
with  the  results  from  previous  studies  conducted  by  the  Sacramento  District. 
These  contour  plots  w'ere  generated  using  nodal  values  of  the  total  normal 
stresses  (that  is,  including  static  stresses).  To  highlight  the  spatial  distribution  of 
tensile  demands,  the  contours  were  generated  using  a  minimum  cutoff  value  of 
150  psi.  The  peak  values  of  maximum  normal  vertical  stresses  across  the  section 
are  696  psi  (San  Fernando  [+H-V]),  698  psi  (Imperial  Valley  [-H-V]),  and 
693  psi  (Chalfent  [+H+V]). 

Figure  31  show's  the  distribution  of  maximum  normal  vertical  stresses  along 
the  upstream  face  of  the  monolith  for  the  San  Fernando  earthquake  (case  [+H-V]). 
To  assess  the  influence  of  the  vertical  component  of  the  input  ground  motion,  two 
different  cases  are  considered:  a)  horizontal  component  only,  and  b)  both  horizontal 
and  vertical  components.  In  addition,  total  (Sy,  including  static  effects)  and 
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Figure  21 .  Maximum  principal  stress  -  monolith  21  (upstream)  -  San  Fernando 
earthquake  [+H-V] 


Figure  22.  Maximum  principal  stress  -  monolith  21  (heel)  -  Imperial  Valley  earthquake 
[-H-V] 
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5  23.  Maximum  principal  stress  -  monolith  21  (downstream)  -  Imperial  Valley 
earthquake  [-H-V] 


Figure  24.  Maximum  principal  stress  -  monolith  21  (upstream)  -  Imperial  Valley 
earthquake  [-H-V] 
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Figure  25.  Maximum  principal  stress  -  monolith  21  (heel)  -  Chalfant  earthquake 
[+H+V] 


Figure  26.  Maximum  principal  stress  -  monolith  21  (downstream)  -  Chalfant 
earthquake  [+H+V] 
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Figure  27.  Maximum  principal  stress  -  monolith  21  (upstream)  -  Chalfant  earthquake 
[+H+V] 


Figure  28.  Normal  vertical  stress  Sy-  monolith  21  - 
San  Fernando  earthquake  [+H-V] 
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Figure  29.  Normal  vertical  stress  Sy  -  monolith  21  - 
Imperial  Valley  earthquake  [-H-V] 


Figure  30.  Normal  vertical  stress  Sy  -  monolith  21  - 
Chalfant  earthquake  [+H+V] 
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dynamic  (Sy ,  excluding  static  effects)  normal  stresses  are  shown  in  the  figure 
(blue  and  red  lines,  respectively).  For  the  horizontal  excitation  case,  the  highest 
stresses  occurred  at  about  72  percent  of  the  height  (Sy  =  574  psi,  =  632  psi). 

When  both  ground  motion  components  are  considered,  the  highest  stresses 
occurred  at  the  base  (Sy  =  696  psi,  Sy  =  761  psi).  Figure  32  shows  the  variation  of 

maximum  normal  stresses  along  the  upstream  face  for  the  Imperial  Valley 
earthquake  (case  [-H-V]).  For  the  case  of  horizontal  excitation  only,  the  highest 
stresses  occurred  at  the  base  (Sy  =  519  psi,  S^y  =  584  psi).  However,  in  the  case  of 

both  horizontal  and  vertical  components,  the  highest  total  stress  occurred  at 
72  percent  of  the  height  (Sy  =  698  psi),  but  the  highest  dynamic  stress  was  found 
at  the  base  (Sy  =  782  psi).  Finally,  similar  plots  depicting  the  maximum  normal 

stresses  along  the  upstream  face  for  the  case  of  the  Chalfant  earthquake  (case 
[+H+V])  are  shown  in  Figure  33.  The  highest  normal  stresses  occurred  at  the 
base  for  both  cases  “a”  (Sy  =  640  psi,  S^  =  706  psi)  and  “b”  (Sy  =  693  psi,  S" 

=  758  psi).  ’ 

Results  -  monolith  14 

As  described  before,  three  different  finite-element  models  were  initially 
considered  for  the  middle  cross  section  of  monolith  14.  The  models  attempted  to 
account  for  the  presence  of  the  spillway  pier  and  the  openings  corresponding  to 
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re  32.  Distribution  of  maximum  total  (blue)  and  dynamic  (red)  normal  vertical  stresses  along  the 
upstream  face  of  monolith  21  -  Imperial  Valley  earthquake  -  cases  [-H]  and  [-H-V] 


Figure  33,  Distribution  of  maximum  total  (blue)  and  dynamic  (red)  normal  vertical  stresses  along  the 
upstream  face  of  monolith  21  -  Chalfant  earthquake  -  cases  [+H]  and  [+H+V] 
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the  chambers  and  conduits.  The  first  model  (Model  1)  accounts  for  the  presence 
of  the  spillway  pier  by  using  equivalent  material  properties;  the  second  model 
(Model  2)  represents  the  spillway  pier,  both  chambers,  and  both  conduits  by 
means  of  equivalent  material  properties;  and  the  third  model  (Model  3)  uses 
equivalent  material  properties  to  represent  the  spillway  pier  and  the  conduits  but 
it  incorporates  the  voids  corresponding  to  the  chambers.  The  motivation  was  to 
capture — -at  least  globally — ^the  effects  of  those  features  on  the  overall  dynamic 
characteristics  of  the  system.  The  proper  evaluation  of  the  state  of  stress  and 
deformation  near  those  features  obviously  requires  three-dimensional  modeling. 
However,  it  was  expected  that  these  approximate  models  would  be  able  to  gener¬ 
ate  reasonable  estimates  of  the  overall  level  of  stress  demands  across  the  section 
of  the  monolith. 

Tables  1 1  and  12  show  the  first  ten  natural  frequencies  corresponding  to  the 
three  models  of  monolith  14  for  empty  reservoir  conditions.  Table  1 1 
corresponds  to  the  rigid  foundation  case  {E^-ca),  and  Table  12  corresponds  to 

the  flexible  foundation  case  (Ej-  =  7.9  •  10®  psi).  Depending  on  the  model,  the 

incorporation  of  foundation  flexibility  causes  an  increase  of  about  35  to 
approximately  39  percent  in  the  fundamental  period  of  the  dam-foundation 
system  with  respect  to  the  dam  on  rigid  base.  In  general,  foundation  flexibility 
increases  all  the  natural  periods,  with  a  higher  influence  on  the  period 
corresponding  to  the  second  vibration  mode. 


Table  11 

Natural  Frequencies  -  Monolith  14  with  Rigid  Foundation  (Empty 
Reservoir) 

Mode 

Rigid  Foundation  (1) 

Rigid  Foundation  (2) 

Rigid  Foundation  (3) 

Frequency 

(Hz) 

Period 

(sec) 

Frequency 

(Hz) 

Period 

(sec) 

Frequency 

(Hz) 

Period 

(sec) 

1 

6.26 

0.160 

6.16 

0.162 

6.06 

0.165 

2 

14.17 

0.071 

14.09 

0.071 

13.68 

0.073 

3 

15.07 

0.066 

14.97 

0.067 

14.80 

0.068 

4 

22.90 

0.044 

22.91 

0.044 

22.20 

0.045 

5 

31.62 

0.032 

31.60 

0.032 

30.13 

0.033 

6 

33.69 

0.030 

33.67 

0.030 

32.34 

0.031 

7 

37.58 

0.027 

37.57 

0.027 

34.92 

0.029 

8 

39.81 

0.025 

39.79 

0.025 

36.58 

0.027 

9 

43.01 

0.023 

43.01 

0.023 

40.16 

0.025 

10 

45.98 

0.022 

45.75 

0.022 

40.98 

0.024 

Based  on  the  objectives  of  the  analyses.  Model  2  was  adopted  as  the  most 
appropriate  representation  of  the  monolith  section.  This  model  incorporates  the 
appropriate  modifications  in  the  material  properties  to  reflect  the  presence  of  the 
spillway  pier,  the  chambers,  and  the  conduits,  while  assuming  an  idealized  2D 
section  of  uniform  thickness.  Although  in  a  simplified  manner,  this  model 
accounts  for  the  chambers  and  the  conduits,  and  therefore  it  is  likely  to  generate  a 
better  approximation  of  the  actual  dynamic  response  characteristics  of  the  struc¬ 
ture  than  Model  1  (which  considers  only  the  influence  of  the  spillway  pier).  On 
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the  other  hand,  Model  3  exaggerates  the  influence  of  the  chamber  voids,  whose 
width  (in  the  direction  normal  to  the  plane  of  the  section)  becomes  then  equal  to 
the  monolith  thickness  in  a  2D  representation  that  assumes  a  section  of  uniform 
thickness.  The  locations  of  the  first  five  natural  periods  of  Model  2  are  shown  in 
Figures  34  and  35  with  respect  to  the  MCE  5-percent-damped  horizontal  and  ver¬ 
tical  response  spectra.  The  results  in  Figures  34  and  35  correspond  to  rigid  and 
flexible  foundations,  respectively,  and  they  were  computed  for  empty  reservoir 
conditions. 


0.01 

0.10 

1.00 

Period  [sec] 

Figure  34.  Natural  periods  and  MCE  spectra  -  monolith  14  with  rigid 
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Figure  35.  Natural  periods  and  MCE  spectra  -  monolith  14  with  flexible  foundation 


Figure  36  shows  the  distribution  of  static  normal  vertical  stress  (Sy)  along  the 
dam-foundation  interface,  considering  only  self-weight  and  hydrostatic  loading. 
The  results,  which  correspond  to  Model  2,  are  presented  for  both  rigid  and  flexi¬ 
ble  foundation  conditions.  The  figure  shows  that  for  the  assumed  reservoir  level 
(gross  pool  elevation  at  466  ft),  the  static  normal  stress  at  the  heel  of  the  dam 
reaches  values  of  1 15  psi  and  -48  psi,  for  rigid  and  flexible  base  conditions, 
respectively.  The  distribution  of  the  normal  vertical  stress  Sy  across  the  section  is 
depicted  by  the  contour  plots  shown  in  Figures  37  and  38,  for  rigid  and  flexible 
foundation  cases,  respectively. 

Tables  13  through  15  indicate  the  peak  values  of  the  maximum  principal 
stresses  in  the  section  computed  for  the  San  Fernando,  Imperial  Valley,  and 
Chalfant  earthquakes,  respectively.  The  principal  stress  values  reported  in  these 
tables  correspond  to  element  (centroidal)  stresses,  and  they  include  the  static  con¬ 
tributions.  For  all  three  earthquakes,  the  highest  tensile  stress  occurred  at  the  heel 
of  the  section.  For  the  San  Fernando  earthquake,  the  highest  value  of  maximum 
principal  stress  reached  897  psi,  which  occurred  for  the  [+H-V]  component 
combination.  The  results  corresponding  to  the  Imperial  Valley  and  Chalfant 
earthquakes  indicated  values  of  maximum  principal  stresses  up  to  790  psi  (case 
[-H-V])  and  746  psi  (case  [+H+V]),  respectively.  The  three  critical  component 
combinations  coincide  with  those  identified  for  the  seismic  demand  at  the  base  of 
monolith  21. 

Solely  examining  the  peak  stress  responses  at  selected  locations  does  not  pro¬ 
vide  sufficient  information  to  judge  the  comparative  severity  of  different  ground 
motions,  and  a  thorough  evaluation  of  the  seismic  performance  of  a  concrete  dam 
must  take  into  account  not  only  the  magnitude  of  the  stress  responses  but  also 
their  time-varying  characteristics.  The  time  histories  of  the  critical  maximum 
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Table  15 

Peak  Values  of  Maximum  Principal  Stresses  -  Monolith  14  - 
Chalfant  Earthquake 

Case 

X 

(ft) 

CTmmt  Foundation 
(psi) 

+H 

6.01 

6.07 

9.7 

659 

-H 

6.01 

6.07 

9.9 

693 

+H+V 

6.01 

6.07 

6.0 

746 

+H-V 

6.01 

6,07 

9.7 

682 

-H+V 

6.01 

6.07 

11.6 

689 

-H-V 

|^|^^|[|| 

6.07 

9.9 

739 

principal  stresses  are  presented  in  Figures  39  through  41  for  the  three 
earthquakes.  As  mentioned  before,  the  San  Fernando  case  shows  the  most 
significant  tensile  peak  (897  psi),  but  high  values  appear  as  isolated  tensile 
excursions  (for  example,  the  maximum  principal  stress  history  exceeds  700  psi 
only  twice). 


Figure  39.  Maximum  principal  stress  -  monolith  14  (heel)  -  San  Fernando  earthquake 
[+H-V] 


It  is  also  important  to  evaluate  the  spatial  distribution  of  the  tensile  demands 
across  the  section.  Contour  plots  of  the  peak  values  of  normal  vertical  stresses  (Sy) 
are  shown  in  Figures  42  through  44  for  the  critical  component  combinations 
corresponding  to  the  San  Fernando,  Imperial  Valley,  and  Chalfant  earthquakes, 
respectively.  This  particular  response  quantity  was  selected  to  facilitate  the  com¬ 
parison  with  the  results  from  the  previous  response-spectrum-based  seismic  analy¬ 
ses  performed  by  the  Sacramento  District.  The  stress  contoure  were  generated  using 
nodal  normal  stress  values  that  included  the  effects  of  the  static  loads.  Only 
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Figure  40.  Maximum  principal  stress  -  monolith  14  (heel)  -  Imperial  Valley  earthquake 
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Figure  42.  Normal  vertical  stress  Sy  -  monolith  14  - 
San  Fernando  earthquake  [+H-V] 


Figure  43.  Normal  vertical  stress  Sy  -  monolith  14  - 
Imperial  Valley  earthquake  [-H-V] 
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Figure  44.  Normal  vertical  stress  Sy  -  monolith  14  - 
Chalfant  earthquake  [+H+V] 

tensile  demands  exceeding  150  psi  are  presented  in  the  figures.  The  normal  verti¬ 
cal  stresses  across  the  section  reach  peak  positive  values  of  789,  693,  and  650  psi 
for  the  San  Fernando  (case  [+H-V]),  Imperial  Valley  (case  [-H-V]),  and  Chalfant 
(case  [+H+V])  earthquakes,  respectively.  In  all  cases,  the  peak  values  occurred  at 
the  heel  of  the  section. 

The  influence  of  the  vertical  component  of  the  input  ground  motion  on  the 
stress  response  is  investigated  in  Figures  45  through  47.  They  show  the  peak  val¬ 
ues  of  the  maximum  vertical  stresses  along  the  upstream  face  of  the  dam  section. 
Each  figure  shows  the  stresses  due  to  horizontal  excitation  (case  a)  and  both 
ground  motion  components  (case  b).  Also,  the  influence  of  the  static  stresses  on 
the  total  response  is  evaluated  by  displaying  not  only  the  total  vertical  stresses 
(Sy),  but  also  the  dynamic  vertical  stresses  (S^ ,  excluding  static  effects).  As  previ¬ 
ously  indicated  by  the  contour  plots,  the  highest  tensile  demands  take  place  at  the 
base  for  all  cases.  Figure  45  shows  the  profile  of  vertical  stresses  for  the  San  Fer¬ 
nando  earthquake,  with  peak  values  of  S^,  =  789  psi  and  S^  =  838  psi  for  the  case 

in  which  both  ground  motion  components  act  simultaneously.  Notice  that  if  only 
the  horizontal  component  of  the  earthquake  motion  is  considered,  the  correspond¬ 
ing  peak  values  are  approximately  80  percent  (S”  =  616  psi  and  S™  =  665  psi)  of 
those  obtained  for  both  excitation  components.  Figure  46  corresponds  to  the 
Imperial  Valley  earthquake,  and  it  exhibits  peak  stresses  of  Sy  =  693  psi  and  Sy 
=  742  psi  (both  ground  motion  components).  The  horizontal  excitation  case 
exhibits  peak  values  that  are  again  80  percent  (Sy  =  555  psi  and  S™  =  603  psi)  of 
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45.  Distribution  of  maximum  total  (blue)  and  dynamic  (red)  normal  vertical  stresses  along 
the  upstream  face  of  monolith  14  -  San  Fernando  earthquake  -  cases  [+H]  and 
[+H-\/] 


Figure  46.  Distribution  of  maximum  total  (blue)  and  dynamic  (red)  normal  vertical  stresses  along 
the  upstream  face  of  monolith  14  -  Imperial  Valley  earthquake  -  cases  [-H]  and  [-H-V] 
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Figure  47.  Distribution  of  maximum  total  (blue)  and  dynamic  (red)  normal  vertical  stresses  along 
the  upstream  face  of  monolith  14  -  Chalfant  earthquake  -  cases  [+H]  and  [+H+V] 

those  corresponding  to  both  components  considered  simultaneously.  Figure  47 
shows  the  results  corresponding  to  the  Chalfant  earthquake,  with  peak  values  of 
Sy  =  650  psi  and  Sy  =  699  psi  (both  ground  motion  components).  For  this 
particular  case,  the  peak  stresses  corresponding  to  the  horizontal  excitation  acting 
alone  are  slightly  less  than  90  percent  (S”  =  613  psi  and  S“  =  565  psi)  of  those 
computed  considering  both  earthquake  components. 


Comparison  with  simplified  response  spectrum  analyses 

A  main  objective  of  this  study  is  to  compare  the  results  obtained  by  time  his¬ 
tory  analysis  with  the  maximum  response  estimates  previously  determined  by  SPK 
using  a  simplified  response  spectrum  approach.  Figures  48  and  49  show  the  vari¬ 
ation  of  dynamic  vertical  stresses  along  the  upstream  faces  of  monolith  21  and  14, 
respectively.  The  black  curve  in  those  figures  represents  the  maximum  dynamic 
stress  profile  determined  by  the  simplified  response  spectrum  approach  consider¬ 
ing  both  components  of  the  excitation.  These  results  correspond  to  the  middle 
section  of  monolith  21  (Figure  48)  and  to  a  sectional  cut  at  a  distance  of  18.75  ft 
fi'om  the  center  of  monolith  14  (Figure  49).  The  red  curve  represents  the  peak  val¬ 
ues  of  the  dynamic  stresses  obtained  by  time  history  analysis.  The  time  history 
results  for  monolith  21  (Figure  48)  were  obtained  considering  both  components  of 
the  Imperial  Valley  earthquake  (combination  [-H-V]).  The  time  history  results  for 
monolith  14  (Figure  49)  were  obtained  considering  both  components  of  the  San 
Fernando  earthquake  (combination  [-t-H-V]).  The  stress  results  for  monolith  21 
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Figure  48.  Comparison  of  maximum  values  of  dynamic 
normal  vertical  stresses  along  the  upstream 
face  of  monolith  21  obtained  by  simplified 
response  spectrum  analysis  (black)  and 
time-history  analysis  (red) 
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Figure  49.  Comparison  of  maximum  values  of  dynamic 
normal  vertical  stresses  along  the  upstream 
face  of  monolith  14  obtained  by  simplified 
response  spectrum  analysis  (black)  and 
time-history  analysis  (red) 
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show  good  agreement  along  the  height,  and  the  stress  values  obtained  by  time 
history  analysis  are  slightly  higher  than  those  predicted  by  the  simplified 
response  spectrum  approach.  The  results  for  monolith  14,  however,  show  larger 
differences  between  the  maximum  values  predicted  by  the  two  approaches.  For 
this  case,  the  peak  stress  values  obtained  by  time  history  analysis  are  higher  than 
those  calculated  using  the  simplified  response  spectrum  approach.  It  is  important 
to  consider  the  fundamental  difference  between  the  two  underlying  finite  element 
models.  The  time  history  analysis  is  based  on  a  2D  model  that  incorporates  the 
infiuence  of  the  spillway  pier  and  the  chamber  and  conduit  voids  in  an 
approximate  manner  (using  equivalent  material  properties),  whereas  the 
simplified  response  spectrum  analysis  is  based  on  a  3D  representation  of  the 
geometry.  This  is  very  relevant  because  the  time  history  results  directly  reflect 
the  influence  of  the  equivalent  representation  of  the  spillway  pier  at  the  top  of  the 
model,  while  the  results  from  the  simplified  response  spectrum  approach 
correspond  to  a  sectional  cut  at  18.75  ft  from  the  center  of  the  monolith,  not 
directly  below  the  spillway  pier. 

Figure  50  shows  the  variation  of  maximum  dynamic  vertical  stresses  along 
the  width  of  monolith  14  for  two  different  elevations.  The  two  elevations  are 
209.4  and  277.4  ft,  and  they  correspond  to  the  upper  and  lower  conduit  soffits, 
respectively.  The  black  curves  represent  the  stress  distribution  predicted  by  the 
simplified  response  spectrum  approach  based  on  a  3D  model  of  the  monolith. 

The  stress  increase  induced  by  the  presence  of  the  conduit  openings  is  quite 
noticeable.  The  peak  values  of  dynamic  vertical  stresses  at  those  elevations 
obtained  by  time  history  analysis  are  also  shown  in  the  figure,  and  they  were 
obtained  considering  both  components  of  the  San  Fernando  earthquake 
(combination  [+H-V]).  These  results  were  determined  based  on  a  2D  model  and 
they  are  represented  as  constant  values  across  the  width  of  the  section.  In  both 
cases,  the  peak  values  obtained  by  time  history  analysis  are  bracketed  by  the 
values  determined  by  the  simplified  response  spectrum  approach  across  the 
section  width. 
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Figure  50.  Maximum  dynamic  normal  vertical  stresses  along  the  width  of  the  upstream  face  of 
monolith  14  obtained  by  simplified  response  spectrum  analysis  (black)  and  time- 
history  analysis  (red) 
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4  Conclusions 


This  report  presents  a  series  of  finite-element  based  dynamic  stress  analyses 
performed  on  monoliths  14  and  21  of  Folsom  Dam.  The  linear  time-history 
response  analyses  were  performed  using  the  computer  program  EAGD-84  as  the 
analysis  platform.  The  corresponding  formulation  accounts  for  dam-foundation 
interaction  and  water  compressibility  effects,  and  this  is  entirely  consistent  with 
the  previous  seismic  studies  performed  by  the  Sacramento  District  using 
Chopra’s  simplified  analysis  procedure.  The  middle  cross  sections  of  the 
monoliths  were  idealized  by  two-dimensional  finite-element  models  consisting  of 
four-node  quadrilateral  elements. 

The  time-history  response  analyses  were  performed  using  three  sets  of  input 
motions  spectrally  matched  to  the  adopted  MCE  horizontal  and  vertical  spectra. 

It  is  interesting  to  note  that  in  general  these  three  input  motions  generated 
quantitatively  and  qualitatively  similar  results  because  of  the  excellent  spectral 
matching  of  the  input  time  histories  used  in  the  linear  analyses.  In  all  cases,  the 
critical  loading  scenario  corresponded  to  a  case  in  which  the  base  excitation 
consisted  of  a  combination  of  both  horizontal  and  vertical  input  motion 
components. 

These  two-dimensional  analyses  identified  three  critical  regions  in  the 
nonoverflow  monolith  analyzed  (monolith  21):  the  upstream  dam-foimdation 
contact  comer  (heel),  as  well  as  the  regions  along  both  upstream  and  downstream 
faces  located  at  about  72  percent  of  the  height  of  the  monolith.  The  reported 
principal  stresses  reached  a  peak  value  of  about  760  psi  for  the  San  Fernando 
earthquake.  It  must  be  mentioned  that  this  value  is  representative  of  a  finite  area 
in  the  critical  region  that  is  characterized  by  the  smallest  contributing  element.  In 
the  case  of  a  comer  point,  such  as  the  heel,  there  could  be  significant  stress 
concentrations  induced  by  the  geometry  of  the  section  and  any  imposed  boundary 
conditions.  In  this  case,  the  determination  of  the  peak  stress  value  will  be 
drastically  affected  by  the  local  density  of  the  finite-element  mesh.  Therefore,  it 
is  of  paramount  importance  to  evaluate  not  only  the  peak  value  but  also  the 
corresponding  local  stress  gradient.  Careful  engineering  judgment  should  be 
considered  when  evaluating  these  types  of  stress  concentration  effects.  In  all 
cases,  the  high  tensile  stress  values  w'ere  confined  to  limited  regions  in  the 
neighborhood  of  the  critical  locations.  The  corresponding  stress  time  histories 
indicated  that  critical  tensile  demands  occurred  at  very  isolated  instances  along 
the  duration  of  the  earthquake. 
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The  analyses  perfomed  on  the  two-dimensional  overflow  section  (mono¬ 
lith  14)  identified  the  heel  as  the  critical  region,  where  peak  principal  stresses 
reached  a  value  of  about  897  psi,  also  for  the  San  Fernando  earthquake.  No  other 
critical  areas  were  identified  based  on  the  tj^re  of  finite-element  model  used  for 
these  studies.  Additional  investigations  of  stress  concentration  effects  induced  in 
the  perimeter  of  the  chambers  and  the  conduits  w^ould  require  the  three- 
dimensional  modeling  of  the  monolith  geometry.  As  previously  indicated,  the 
tensile  demands  must  be  evaluated  not  only  as  isolated  peak  values  but  also  con¬ 
sidering  their  temporal  variation  and  the  corresponding  spatial  gradient.  The 
stress  time  history?  at  the  critical  location  indicated  a  few?  isolated  excursions 
above  an  assumed  strength  threshold  (700  psi),  and  in  all  cases  the  areas  of 
excessive  tensile  demands  were  limited  to  the  region  immediately  near  the  heel. 

All  compressive  stresses  predicted  by  these  analyses  are  w?ell  below?  the 
corresponding  strength  limit.  The  levels  of  tensile  demands  predicted  by  this 
study  indicate  that  localized  damage  in  the  form  of  tensile  cracking  is  likely  to 
take  place  in  limited  areas  of  monoliths  14  and  21  when  subjected  to  extreme 
seismic  loading,  but  this  is  still  acceptable  seismic  performance  under  MCE 
conditions. 
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